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OPERATIONAL FREQUENCY STABILITY OF RUBIDIUM AND CESIUM 

FREQUENCY STANDARDS 

John E. Lavery 
Goddard Space Flight Center 


INTRODUCTION 

In the course of testing various rubidium and cesium frequency standards under operational con- 
ditions for use in NASA tracking stations, about 55 unit-years of relative frequency measurements for 
averaging times from 10 to 10 7 s have been accumulated at Goddard Space Flight Center (GSFC). 
Statistics on the behavior of rubidium and cesium standards under controlled laboratory conditions 
have been published by many institutions (see, for example, ref. 1 ), but it was not known to what 
extent the lesser controlled environments of NASA tracking stations affected the performance of the 
standards. The purpose of this report is to present estimates of the frequency stability of rubidium 
and cesium frequency standards under operational conditions based on the data accumulated at GSFC. 


Table 1. -Atomic Frequency Standards Used 
in Experiments 


Serial no. or 
designation 

Manufacturer 

Rb 107 

Varian Associates 

Rb 136 

Varian Associates 

Rb 138 

Varian Associates 

Cs 110 

Hewlett-Packard Co. 

Cs 136 

Hewlett-Packard Co. 

Cs 137 

Hewlett-Packard Co. 

Cs 138 

Hewlett-Packard Co. 

Cs 139 

Hewlett-Packard Co. 

Cs 152 

Hewlett-Packard Co. 

Cs 182 

Hewlett-Packard Co. 

Cs 185 

Hewlett-Packard Co. 

Cs 186 

Hewlett-Packard Co.' 

HM: 


H-10no.2 

Varian Associates 

NX-1 

e) . 


a An experimental hydrogen maser developed at GSFC. 
See ref. 2. 


DATA DESCRIPTION 

The three rubidium gas cells (designated Rb) 
and nine cesium beam frequency standards (Cs) 
on which the measurements were made, as well 
as the two hydrogen masers (HM) used as refer- 
ences for many of the tests, are listed in table 1 
along with their serial numbers or designations 
and their manufacturers. During the tests the 
standards were kept in a laboratory at GSFC. 
Except for the shielding built into the standards 
themselves, there was no special control of the 
ambient magnetic, electric, vibration, and tem- 
perature conditions. The ambient magnetic and 
electric conditions were typically noisy. The 
standards were driven by ac power and were in 
no way isolated by transformers. Vibration from 
nearby air conditioning equipment and from 
trucks at a nearby loading platform was not 
shielded in any way. The ambient temperature 
was typically between 298 and 303 K. There 
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were, however, several brief excursions to temperatures as low as 29 1 K and as high as 3 13 K due to 
equipment problems. These conditions are less controlled than those in the NASA tracking stations. 
Hence the stabilities of the standards when operating in the tracking stations should be at least as good 
as the stabilities calculated in this paper. 

The measurements made on the standards consisted of average relative frequency measurements 
for varying averaging times. In some of the data sets, average relative frequency measurements were 
missing or were bad because of ac power failure or recorder failure. All such points were a posteriori 
linearly interpolated from the nearest earlier (in epoch time) good average relative frequency measure- 
ment and the nearest later (in epoch time) good average relative frequency measurement. 

The total number of measurements made for all types of data used in this report are given in 
table 2. Data sets are said to be of the same “type” when the following parameters are the same for 


Table 2.— Average Relative Frequency Data Sets 


Type of data 

Number of 
data sets m 

Number of measurements 

Test 

unit 

Reference 

unit 

Averaging 

time 

?0’ s 

Dead 
time 
d, s 

Total a 

Interpolated 

Rb 

Rb 

3 600 

0.0 

2 

3 090 

0 

Rb 

Cs (10-s TC) 

10 

2.3 

1 

1 076 

18 

Rb 

Cs (10-s TC) 

100 

2.2 

1 

538 

1 

Rb 

Cs (10-s TC) 

1 000 

2.7 

1 

223 

0 

Rb 

HM 

10 

2.3 

13 

8 473 

67 

Rb 

HM 

100 

2.2 

10 

6 405 

16 

Rb 

HM 

1 000 

2.7 

9 

5 126 

15 

Rb 

HM 

3 600 

.0 

7 

13 320 

308 

Cs 

Cs 

3 600 

.0 

3 

8 851 

263 

Cs (10-s TC) 

HM 

10 

.2 

8 

4 841 

0 

Cs (10-s TC) 

HM 

100 

.2 

8 

4 871 

11 

Cs (10-s TC) 

HM 

1 000 

.2 

8 

4 787 

25 

Cs (60-s TC) 

HM 

10 

.2 

8 

4 634 

0 

Cs (60-s TC) 

HM 

10 

2.3 

3 

1 904 

2 

Cs (60-s TC) 

HM 

100 

.2 

8 

4 706 

0 

Cs (60-s TC) 

HM 

100 

2.2 

3 

2 496 

0 

Cs (60-s TC) 

HM 

1 000 

.2 

8 

4 804 

3 

Cs (60-s TC) 

HM 

1 000 

2.7 

1 

692 

0 

Cs 

HM 

3 600 

.0 

13 

37 404 

1391 

Cs 

HM 

604 800 

.0 

i 

1 

88 

6 


TC = time constant. 

a Total number of measurements for all m data sets, including the interpolated measurements. 
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each set: test unit, 1 reference unit, duration or averaging time r 0 of each average relative frequency 
measurement, and dead time d between successive measurements (that is, the time during which no 
measurement was taken). The servo time constants are indicated only for the cesium standards and 
only when r Q < 1000 s. The difference in effect of a 10- and 60-s time constant for r 0 > 3600 s can 
be neglected because the time constants in such cases are too small with respect to t q to have an 
appreciable effect. The rubidium standards tested all have a fixed servo time constant which is on the 
order of 1 ms. 

Neither temperature effects nor long-term frequency drift was removed from the data before 
analysis because the object of the tests was to measure the stability of the frequency standards under 
operational conditions, where both temperature fluctuations and long-term frequency drift are present. 


STATISTICAL ANALYSIS 


Let there be given a set of m identical test frequency standards and a set of m identical reference 
frequency standards. Let <t> n (f), 1 < n < m, denote the instantaneous fluctuations (measured in time 
units) of the epoch time output of the nth test standard compared to the epoch time output of the nth 
reference standard. Let y n (?) be the instantaneous (fractional) frequency fluctuation of the nth test 
standard compared with the nth reference standard; i.e., 




d<t> n {t) 

dt 


( 1 ) 


Let y„(t) be the average relative (fractional) frequency fluctuation of thenth test standard compared 
with the nth reference standard: 


1 rt+r <P n (t + t )- <j> n (t) 

y n (0= - I y„ (0 dt = 

T Jt T 


( 2 ) 


The constant r is called the averaging time of y(t). The Allan standard deviation a( 2, T, r) of the fre- 
quency fluctuations of the set of test standards compared with the set of reference standards is defined 
to be (ref. 3) 


o(2, T, t) = 



<var [y n (t+T)~ y n (t)]) 


(3) 


where the symbol < > denotes infinite epoch time average. The analysis of all data listed in table 2 con- 
sisted in the calculation of an estimate, which is denoted by s(2, T, r), of a(2, T, r) in the following 
manner. 

Taking any type of data from table 2, let the number of average relative frequency measurements 
in the nth data set, 1 <n<m, bem n . Denote this nth set of average relative frequency measurements 


^though there are sometimes significant differences in the frequency stabilities of various rubidium standards, the three 
rubidium standards listed in table 1 all had mutually close stabilities. For this reason, these rubidium standards will be considered to be 
identical. Because the nine cesium standards listed in table 1 all had mutually close stabilities, they too will be considered to be identical. 


3 



by |>Vi (*)}/=" ■ For i = 1,2 , ,m n - 1, denote the variance of the two average relative frequency 

measurements and y n {i + 1) by v n (i): 


[y„0'+ l)- y„0')l 2 

V n (i)= 


( 4 ) 


The square root of the average over both / ( 1 < i < m n - 1) and n ( 1 < n < m) of these v n (i) is the 
desired estimate of a( 2, r Q + d, r 0 ): 


s( 2, r n + d, r n ) = 


/ rn 

m n~ 1 


I 

I 


/ n = 1 

i= 1 

f 


- 1) 

f n - 

i 



(5) 


From the original data sets {y n \ , 1 < n < m, new data sets with averaging time = 2t q and 
dead time d (assumed small with respect to r 0 ) can be approximated by defining 


y„0';i) = 


y n 0'+ l) + y n ( 0 


( 6 ) 


i = 1 , 2, . . . ,m n - 1 and n= 1,2, ... ,m. Denote the variance of y n (i; 1) and y n ii +2; 1) by v n (i; 1): 

[y n (/ + 2; 1)- y n (i; l)] 2 


v„(i; 1) = 


(7) 


i= 1,2, ... ,m n - 3 and n= 1,2, ... ,m. Estimate a( 2, t x +d,r x ) by 2 



( 8 ) 


Let k be the exponent of the largest power of 2 contained in any of the m n , 1 < n < m. For 
j= 2,3, . . . , k- 1, the data set \y n (i',j)) "^ _2/ +1 with averaging time t. = 2h Q and dead time d is 
successively calculated from the data set \y„(i,j - 1)}”" _2/ +1 by pairwise averaging: 

y (i + 2 ! ~ 1 ; / - l)+y n (i;j- 1 ) 

y„ (/;/) = " ( 9 ) 


2 Throughout this paper the convention is adopted that whenever a summand, e.g., m n - 3 in ~ 3 ), is less than zero, it is 

treated as zero; and whenever a summation, e.g., v n (i\ 1), has an upper limit that is less than the lower limit, it also is treated as 

zero. 
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i= 1, 2, . . . , m n - 2' + l;n = 1, 2, . . . , m;/fixed. Denote the variance of y n (i;j) andy n (/ + 2 ; ;/) by 




[ y „ 0 ' + 2 / ;/)- (/;/)] 2 

2 

/ = 1 , 2, . . . , m n - 2 /+1 + 1 and n = 1 , 2, . . . , m. Estimate a(2, r ; . + d, 7y) by 


5(2, Tj + d, r ; .) 


An example of this procedure for zero dead time is presented in figure 1 . The quantity v repre- 
sents the variance between the ordinates of the two lines to which the dotted line near v points. 


/ - 

m - 2 /+ 1 + 1 


I 

I 

v n (/;/) 

/"“ 1 

/= 1 

/ 1 

(% - 

1 + 1) 

J 71 = 1 




Figure 1. -Calculation of s(2, x, r). 
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RESULTS 


For each type of data listed in table 2 and for each averaging time t. = 2 !t q ,0</<k - 1 (r Q and 
k change with the type of data), the estimate s(2, r. + d, r ; .) of a(2, t. + d, t.) was calculated. 3 The 
results are presented in table 3 and figure 2 for all data involving a rubidium standard as either the test 
or the reference unit and in table 4 and figure 3 for the cesium versus cesium and cesium versus hydro- 
gen maser data. 

In order to use the data in tables 3 and 4 to estimate the frequency stability of the rubidium and 
cesium standards tested, rather than the relative frequency stability of a comparison of two of these 
standards or of a comparison of one of these standards to a hydrogen maser, the following procedure 
is used. Denote the Allan standard deviations of the test standard versus a hypothetical perfect 
standard, the reference standard versus a hypothetical perfect standard, and the test standard versus 
the reference standard by o T ( 2, t + d, r), o R (2, r + d, t), and o TR (2, r + d, r), respectively. Because 
the variances a\{2, r + d, r) and o R (2, t + d,r) are linear functions (in fact, weighted integrals) of the 
respective power spectral densities of the test and reference standards (ref. 3), and because the power 
spectral density of the comparison of two frequency standards is the sum of the power spectral densi- 
ties of each of the standards, the following relation occurs: 

o^(2, r + d, r) = a^,( 2, r + d, t) + o R (2, r + d, r) (12) 

For comparisons of two identical standards (rubidium standard versus rubidium standard and 
cesium standard versus cesium standard), o R (2, r + d, t) = a T ( 2, r + d, t). Hence, from relation (12), 

®T-R r + d, t) 

M2,r + d,r)= (13) 

a/Z 

For all data for which a hydrogen maser was used as a reference, it is assumed that the instabilities of 
the maser were sufficiently small so as to have 


a r (2, r + d, r) « a TR (2, r + d, r) 

The normalized standard deviation a T ( 2, t, r) can be calculated from a T (2, t + d,r) by the 
relation 


a r (2, t + d,r) 

o t (2, r, r) = 

\/B 2 (r,n) 


(14) 


(15) 


where B 2 (r, n) is a bias function (defined in ref. 4); r = (r + d)/r; and fx, representing the type of noise 
of the standard for the fixed averaging time t and fixed dead time d, is determined from 

o r (2, t + d, t) a W 2 (16) 


3 The analysis was carried out by programs E00016 and E00036 of the GSFC Computer Program Library. Program E00016 is 
for input relative phase data; program E00036 is for input relative frequency data. Although program E00016 reads relative phase data 
as input, its output is the Allan standard deviation of relative frequency s(2,?y + <2, rj) defined in eqs. (5), (8), and (11). These two 
programs are based on a program written by David W. Allan of the National Bureau of Standards, Boulder, Colo. 
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Table 3.— Rubidium Standard Frequency Stability 
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Table 4.— Cesium Standard Frequency Stability 
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16 000 .2 .729 3 200 2.2 1.931 4 838 400 .0 .181 

32 000 .2 .674 6 400 2.2 1.318 9 676 800 .0 .167 

64 000 .2 .481 12 800 2.2 1.003 19 353 600 .0 .114 

128 000 .2 .459 25 600 2.2 .883 











It is of interest to note that for the r and p of the data analyzed in this report, B 2 (r, p) differs from 
unity by less than 0.1 percent and can be ignored. Hence, for the data in this report, 

o t (2, t, t)^o t ( 2, r + d, t) (17) 

Of course, relation (1 7) is an exact equality whenever d = 0. 

Using the estimates s( 2, r + d,r) of a TR (2, t + d,r) from figures 2 and 3 in relations (13) and 
(14) and using relation (17), the standard deviations o T ( 2 , t, r) of the rubidium and cesium standards 
tested can be estimated. These estimates of o T ( 2, r, t) are presented in figure 4 as the “operational 
environment” curves. Also shown in figure 4 are curves taken from references 1 and 5 representing 
the performance of rubidium and cesium standards in a “controlled environment.” By “controlled 
environment” is meant an experimental environment shielded from magnetic, electric, vibration, and 
temperature effects much more than the “operational” environment in which the data presented in 
figures 2 and 3 were taken. 4 The upper curve for rubidium standards under a controlled environment 
in figure 4 is taken from reference 5 and represents the measured performance of Varian rubidium 
standards under controlled conditions. The lower curve for rubidium standards under a controlled 
environment and the curve for cesium standards under a controlled environment in figure 4 are taken 
from reference 1 and represent the measured performance of Hewlett-Packard rubidium and cesium 
standards under controlled conditions. 

CONCLUSIONS 

From figure 4 it is apparent that an operational environment degrades the performance of the 
rubidium standards (by up to one order of magnitude) for frequency averaging times between 10 and 
1 0 3 s and that it degrades the performance of the cesium standards (by up to one order of magnitude), 
for frequency averaging times between 3 X 10 4 and 2 X 10 7 s. For all other averaging times in the 
range covered by the data in figure 4, the stabilities of the standards are not degraded by the opera- 
tional conditions. 
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4 That this was the case for the data presented in reference 1 was verified by private communication with G. M. R. Winkler. That 

the environment in which the data presented in reference 5 was taken was “controlled”, is assumed. 
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